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Abstract 

Starting from integrable cellular automata we present a novel form of Painleve equa- 
tions. These equations are discrete in both the independent variable and the dependent 
one. We show that they capture the essence of the behavior of the Painleve equations 
organizing themselves into a coalescence cascade and possessing special solutions. A 
necessary condition for the integrability of cellular automata is also presented. 



A novel extension of the Painleve equations, well-known for their numerous applica- 
tions in mathematics and physics, will be presented in this paper. Recent progress in 
the domain of integrable discrete systems has led to the derivation of discrete analogs 
of the Painleve equations [1] . They were first identified in a 2D model of quantum grav- 
ity where they appeared as an integrable recursion relation for the calculation of the 
partition function [2]. With hindsight, their first occurence can be traced back to the 
work of Jimbo and Miwa on the spin-spin correlation function of a 2D Ising lattice [3] . 
To date, the full list of discrete Painleve equations has been established [4] and their 
properties (in perfect parallel with those of their continuous counterparts) are being 
actively investigated. Moreover the study of discrete integrable systems has made clear 
the fact that they, rather than the continuous ones, are the fundamental entities. In 
fact, the latter are contained, through appropriate limits, in the discrete ones and thus 
we can establish the hierarchical diagram: 

Discrete systems — > Semi-continuous systems — > Continuous systems 

A recent discovery has made possible to extend this structure in another direction. 
In [5] one of the present authors has proposed a systematic way for the introduction 
of integrable ultra-discrete systems. In previous works on discrete systems, while the 
independent variables were discrete, the dependent variables were assumed to vary con- 
tinuously. The ultra-discrete limit provides a systematic way to discretize the depen- 
dent variable. One can, starting from a given evolution equation, obtain the cellular 
automaton (CA) equivalent. The aim of this paper is two-fold. First, introduce the 
ultra-discrete analogs of the Painleve equations and investigate their properties and, 
second, provide integrablility conditions for cellular-automaton like equations. 

Cellular automata have been the object of an impressive number of studies and 
their behavior is known to be of the utmost richness. The integrability of such sys- 
tems has not been thoroughly studied, since it represents considerable difficulties. An 
occurence of an integrable automaton has been noted in [6] by Pomeau who obtained 
explicitly its constant of motion. Cellular automata representing evolution equations 
have been studied from the point of view of the existence of localized, soliton-like, solu- 
tions. The notion of soliton for CA's was first introduced by Park et al. in [7]. Further 
examples of such CA's with soliton-like structures were given by the Clarkson group [8]. 
Integrable CA were introduced by Bruschi and collaborators [9] who derived Lax pairs 
for their cellular-automaton equations. Bobenko et al. [10] have proposed an interest- 
ing approach to integrable CA's by considering them as the restriction of an integrable 
discrete equation over a finite field. However in many cases the relation to the well- 
known integrable evolution equations was based on circumstancial evidence rather than 
a systematic derivation. The situation has changed recently due to the introduction 
of a method [5] that allows one to convert a given discrete evolution equation to one 
where the dependent variable also takes discrete values. The starting point was the 
CA proposed by one of the authors in collaboration with Satsuma [11]. This simple 
model (essentially a "box and ball" system) was shown later to be the ultra-discrete 



limit of the KdV equation [5]. This was obtained through a limiting procedure (the 
ultra-discrete limit) which allows one to derive a cellular automaton equation starting 
from the appropriate form of a discrete evolution equation. 

As an illustration of the method and a natural introduction to ultra-discrete 
Painleve equations, let us consider the following discrete Toda system: 

u t + 1 -2u t n + u t n - 1 =log(l + 5 2 (e<+i-l))-21og(l + 5 2 (e<-l))+log(l + 5 2 (e<-i-l)) 

(1) 

which is the integrable discretization of the continuous Toda system: 

d 2 r n 



dt 2 



e r ™ +1 - 2e r " + e 7 ^" 1 (2) 



We introduce w through 5 = e~ L / 2e , w l n = eu l n — L and take the limit e — > 0. We use 
the well-known result lim e ^ elog(l + e x / e ) = max(0, x) = (x + \x\)/2 = (x) + . The 
function (x)+ is also known under the name of truncated power function and is equal 
to for x < and x for x > 0. Thus the ultra discrete limit of (1) becomes simply [12]: 

< +1 - 2«4 + K' 1 = - 2K)+ + « _i)+ (3) 

Since the truncated power function of any integer is a non-negative integer, it is clear 
that if w has initially integer values the values will remain integer at all subsequent time 
steps. Thus equation (3) is indeed a cellular automaton equation. 

Let us now restrict ourselves to a simple periodic case with period two i.e. r n+2 = 
r n and similarly w n+ 2 = w n . Calling tq = x and r\ = y we have from (2) the equation 
x = 2e y — 2e x and y = 2e x — 2e y resulting to 'x + y = 0. Thus x + y = [it + v and we 
obtain after some elementary manipulations: 

x = ae^e~ x - 2e x (4) 

Equation (4) is a special form of the Painleve Pm equation. Indeed, putting v = e x -f lt / 2 i 
we find: 

v=— + e ^' 2 (a-2v 2 ) (5) 
v 

The same periodic reduction can be performed on the ultra-discrete Toda equation (3). 
We introduce Wq = X 1 : , w\ = Y l and have, in perfect analogy to the continuous case, 
X t+1 - 2X l + X l ~ l = 2(Y t ) + - 2(X t ) + and Y t+1 - 2Y t + Y*' 1 = 2(X t ) + - 2(Y t )+. 
Again, A 2 (X* + Y l ) = and we can take X 1 + Y l = rat + p (where m, t,p take integer 
values). We find thus that X obeys the ultra-discrete equation: 

X t+1 - 2X l + X 1 ' 1 = 2(mt +p-X t )+- 2(X t )+ (6) 



This is the ultra-discrete analog of the special form of the Painleve Pni equation (5). 
Figure 1 gives a comparison of the evolution under (4) and (6). We remark that the 
dynamics are very similar: the two particles converge towards each other, rebound once 



or twice, get captured and go on oscillating around some equilibrium point. Thus, start- 
ing from a well-known physical problem we have introduced the corresponding cellular 
automaton equation and, restricting it to the simplest periodic lattice, we obtained the 
ultra-discrete form of a Painleve equation. 

In order to construct the ultra-discrete analogs of the Painleve equations (u-P) we 
must start with the discrete form that allows the ultra-discrete limit to be taken. The 
general procedure is to start with an equation for x, introduce X through x = e x ^ e and 
then take appropriately the limit e — > 0. Clearly the substitution x = e x ^ requires x 
to be positive. This is a stringent requirement that limits the exploitable form of the 
d-P's to multiplicative ones. Fortunately many such forms are known for the discrete 
Painleve transcendents [13]: 
d-Pj 

a\ n 1 

Xn+iXn-i = h -~ (7a) 

Xn+lXn-! = CtA" + — (76) 
Xn+xXn-X = a\ U X n + 1 (7c) 



d-P 



d-P 



in 



\"(x n + a\") 

Xn+lXn-l = 7T— 7 (7rf) 

Xn\i- "T X n ) 

X n -\- CtA , . 

X n+ iX n -i = — : — [(e) 

Xn+lXn ~ 1 ~ (l + 7 x n A»)(l + fo„A») U/J 

We remark that in some cases, more than one form exists for a given d-P. The derivation 
of the equivalent ultra-discrete forms is elementary: we take the logarithm of both sides 
of the equation and whenever a sum appears we apply the limit leading to the truncated 
power function. We find thus: 
u-Pi 

X n+l + X n _ x + 2X n = (X n + n + a)+ (8a) 
X n+1 +X n _ 1 +X n = (X n + n + a)+ (86) 
X n+l +X n _ l = (X n + n + a)+ (8c) 



u-Pn 

X n+1 +X n _ 1 =n + (n + a-X n ) + - (X n )+ (8d) 
X n+1 + X n _! - X n = (n + a - X n ) + - (X n + n + 6)+ (8e) 

u-Pni 

X n+ i+X n _i-2X n = (n+a-X n ) + + (n+b-X n ) + -(X n +c+n) + -(X n +d+n)+ (8f) 



These equations describe celular automata provided we restrict the values of the param- 
eters as well as the initial values of the dependent variable to integers. Figure 2 shows a 
comparison of the discrete Pi (7c) with the ultra-discrete Pi (8c). It is remarkable that 
the behavior of the two equations is very similar (provided that we plot the logarithm of 
the variable of the discrete equation, as expected). The forms (8) are not canonical in 
the sense that they can be simplified by a translation of X and a linear transformation 
in n. We will not enter into these details but merely list the canonical forms obtained: 
canonical u-Pi 

X n+1 +X n _i + aX n = (X n + n)+ with a = 0,1,2 (9a) 

canonical u-Pn 

X n+1 +X n - 1 -aX n = a + (n-X n ) + -(n + X n ) + with a = 0, 1 (96) 
canonical u-Pm 

X n+ i+X n _i-2X n = (n + a-X n ) + + (n-a-X n ) + -(X n +b+n) + -(X n -b+n) + (9c) 

Note that equation (6), for m = 2,p = is the subcase a = b = of (9c) after the 
change of variable X 1 = X n + n. 

At this point two questions appear unavoidable. First, is it justified to call these 
equations ultra-discrete Painleve equations? What do they have in common with the 
familiar Painleve equations? One first remark is that the u-P's form a coalescence 
cascade just like their continuous and discrete counterparts [13]. Indeed, starting from 
u-Pn (8d) (respectively (8e)) we can take b — > — oo in which case the nonlinear term 
containing it is always zero. We transform the resulting equation using the identity 
(x)+ = x + (— x)+, then translate X and obtain equation u-Pi in the form (8b), (resp. 
(8c)). We can also recover (8a) from (8d) (and (8b) also from (8c)). Similarly starting 
from u-Pm (8f) we can obtain u-Pn (8d) by taking b — > +oo and c — > +oo such that b—c 
is finite. Next we translate X and through a linear transformation of n we find u-Pn 
(8d). This is not the only property the u-P's share with the continuous and discrete 
Painleve equations as we shall see below. The second question is whether it is possible 
to guess the forms of the u-P's. In other words, what is the (integrability) criterion that 
would single them out among all possible equations? 

In the case of discrete systems the criterion for integrability (equivalent to the 
Painleve property) is the property known as singularity confinement [14]. For cellular 
automata no singularity can exist and thus this criterion is inoperable. The situation is 
analogous to polynomial mappings where no singularity can appear. There, the criterion 
for integrability is based on arguments of growth of the degree of the iterate (or the 
similar notion of complexity introduced by Arnold [15]). Veselov [16] has shown that 
among mappings of the form x n+ i — 2x n + x n -\ = P(x n ) with polynomial P(x), only 
the linear one has non-exponential growth of the degree of the polynomial that results 



from the iteration of the initial conditions. Let us apply such a low-growth criterion to 
a family of ultra-discrete Pi equations of the form: 

X n+1 + aX n + X n _x = (X n + <f>(n))+ (10) 

The three u-Pi obtained from (7) correspond to a = 0, 1, 2. What is the condition for X 
not to grow exponentially towards ±oo? We ask simply that the polynomials r 2 +ar + 1 
and r 2 + (a — l)r + 1 have complex roots (otherwise exponential growth ensues). The 
only integer values of a satisfying this criterion are a = —1,0, 1,2. We remark that 
the three values mentioned above are exactly retrieved plus the value a = — 1. A close 
inspection of this mapping shows that it is also integrable: it is just a form of an ultra- 
discrete Pin, obtained from the discrete system x n+ ix n -i = x n (x n + X n ) which leads 
to (10) with <f>{n) = 0. 

We have applied the low-growth criterion to other cases like u-Pn and u-Pm and 
in every case the results of the growth analysis correspond to the already obtained in- 
tegrable cases. However low-growth is not a sufficiently powerful integrability criterion. 
In particular the inhomogeneous terms ((f) in equation (10)) cannot be fixed by this argu- 
ment. Any <f>(n) would satisfy this requirement. So another criterion must complement 
this first one. 

In the case of (continuous) evolution equations two integrability criteria are often 
used in conjunction: the Painleve property and the existence of multisoliton solutions. 
In the case of Painleve equations the latter are the special solutions that exist for 
particular values of the parameters [17] (except for Pi which is parameter-free). A 
particular class of these solutions (existing also in the discrete case) are the rational ones. 
We shall investigate this property in the case of u-P equations. This will strengthen 
the argument that (9) are indeed Painleve equations and will allow us to fix completely 
their form. For d-Pn the simplest rational solution is a constant. Thus a constant 
solution should exist for u-Pn and indeed for (9b) with a = we find that X = is a 
solution for a = 0. However this solution exists if we replace n in (9b) by any function 
of n. The next solution for u-Pn is a step-function one. Indeed when n is large constant 
positive solution exists with X equal to a/4 while a constant solution with X = a/2 
exists when n is large negative. Thus, for instance for a = 4 a solution for n « is 
X = 2, while for n » a solution is X = 1. The remarkable fact is that that these 
constant "half solutions do really join to form a solution of (9b) with a unique jump 
at n = — 1. It is straightforward to check that this will not be the case in general if the 
non-autonomous part is not linear in n. The general solution of this type becomes now 
clear. For a = Am we have a solution with m jumps from the value X = 2m to X = m. 
The first jump occurs at uq = 1 — 2\m\ and we have successive jumps of —\m\/m at 
n = n + 3/c, k = 0, 1, 2, ... , \m\ — 1. Analogous results can be obtained for the other 
u-Pn corresponding to a = 0. Thus u-Pn has a rich structure of particular solutions. 

On the basis of these results we can conclude that the integrability criterion for 
automata-like Painleve equations appears to be based on low-growth arguments together 



with the existence of a rich class of explicit, globally described, solutions. The need for 
such a two-step process is the price we have to pay in the CA case because of the absence 
of singularities. Can these ideas be applied to other CA evolution equations? Let us 
present the example of ultra-discrete Burgers equation: 

Using the identity (X) + = X + (—X) + we can rewrite (11) as: 

X t + 1 = a(X t n+1 ) + -b(-X t n ) + (12) 

where, from (11), a = b = 1. How can we obtain these values of the parameters based 
on a low-growth argument? If either \a\ > 1 or \b\ > 1 we can find a special direction 
where either X^ +k or X^j^ 1 grows exponentially. Thus the only integer values that 
a and b can take are ±1 and 0. The value leads to uninteresting, one dimensional 
evolution and the case a = b = — 1 is equivalent to the case a = b = 1. Thus the only 
other case that cannot be excluded on low-growth arguments is the case a = 1 and 
b = — 1. Although the discrete equation from which we could have obtained it is not 
integrable the behavior of the automaton-like equation is regular and there is as yet no 
indication as to its nonintegrability. 

It is clear from the results presented in [5] and the ones above that we can pro- 
duce cellular-automaton-like equations starting from integranble nonlinear evolution 
equations. The present work has focused on the ultra-discrete forms of Painleve equa- 
tions of which the first three were given. The forms of the remaining u-P equations 
will also be investigated: probably they will require a two-component description. The 
properties of their solutions as well as the relations between the various u-P's must also 
be studied. These open questions will be addressed in future works. What is important 
at this stage is that we have shown that this new domain of integrable systems is par- 
ticularly rich. While the discrete systems are the fundamental entities and contain all 
the structure, the cellular automata are their bare-bones version capturing the essence 
of the dynamics. This explains the interest that these ultra-discrete systems present 
both from the mathematical and the physical points of view. 
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Figure captions. 

Figure 1. Distance between the two particles as a function of time in the case (a) of the 

continuous Toda potential and (b) its ultra-discrete analog. 
Figure 2. Solution of the discrete Painleve I equation (a) and of its ultra-discrete analog (b). 



